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ABSTRACT 

The chemical composition of the ultra-high-energy (UHE) cosmic rays serves as an important clue 
for their origin. Recent measurements of the elongation rates by the Pierre Auger Observatory hint 
at the possible presence of heavy or intermediate mass nuclei in the UHE cosmic rays. Gamma-ray 
bursts (GRBs) and hypernovae have been suggested as possible sources of the UHE cosmic rays. 
Here we derive the constraints on the physical conditions under which UHE heavy nuclei, if they are 
accelerated in these sources, can survive in their intense photon fields. We find that in the GRB 
external shock and in the hypernova scenarios, UHE nuclei can easily survive photo-disintegration. 
In the GRB internal shock scenario, UHE nuclei can also survive, provided the dissipation radius 
and/or the bulk Lorentz factor of the relativistic outflow are relatively large, or if the low-energy 
self-absorption break in the photon spectrum of the prompt emission occurs above several KeV. In 
internal shocks and in the other scenarios, intermediate-mass UHE nuclei have a higher probability of 
survival against photo-disintegration than UHE heavy nuclei such as Fe. 
Subject headings: gamma rays: bursts — cosmic rays 



1. INTRODUCTION 

Ultra-high-energy (UHE) cosmic rays above ~ 10 18 eV 
are thought to be of extra-galactic origin, since charged 
particles with such high energies cannot be confined by 
the magnetic field of our Galaxy. Galactic cosmic ray ac- 
celerators, such as supernova remnants, are expected to 
reach their maximum energy below ~ 10 18 eV. The tran- 
sition between galactic and extra-galactic cosmic rays is 
therefore believed to happen either at the 'second knee' 
around 10 18 eV, where the chemical composition changes 
significantly as measured by HiRes, followed by a 'dip' in 
the spectrum (Berezinsky et al. 2006), or at the 'ankle' 
around 10 19 eV where the cosmic ray spectrum becomes 
flatter. 

At these ultra- high energies, the chemical composition 
is a subject of debate. It has been speculated that these 
cosmic rays are made of pure protons, up to the highest 
energies. On the other hand, there are also reasons for 
favoring a cosmic ray spectrum dominated by heavy or 
intermediate mass nuclei at the highest energies, since ac- 
cording to the Hillas criterion (Hillas 1984) astrophysical 
sources are able to accelerate particles up to a maximum 
energy proportional to their nuclear charge. Recently, a 
mixed composition scenario has been invoked to fit the 
UHECR spectrum above - 10 185 eV (Allard et al. 2005). 
By studying the details of the development of the result- 
ing air showers, one can in principle infer the species of 
the primary UHECRs, since at a given energy, showers 
initiated by heavy nuclei develop higher up in the at- 
mosphere than proton-induced showers. Both AGASA 
and HiRes data favor a dominance of light hadrons, con- 
sistent with being pure protons, in the composition of 
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UHECRs above 10 19 eV (Hughes et al. 2007), which 
is consistent with models assuming UHECRs above 10 18 
eV are due to extra-galactic protons. On the other hand, 
a recently report of an analysis of the muon component 
of UHE air showers in the Yakutsk data showed a heavy 
nuclei component in the UHECR spectrum (Glushkov et 
al. 2007). 

One of the aims of Pierre Auger Observatory is to 
study the composition of UHECRs. It provides more pre- 
cise measurements of the depth of the UHECR-initiated 
shower maximum Y max at extremely high energies (> 
10 19 eV), with uncertainties which are about a factor 4 
smaller than those of the best measurement from the 
HiRes experiment (Unger et al. 2007). The elonga- 
tion rate data presented by the Pierre Auger Observa- 
tory team is satisfactorily represented by a fit containing 
a break point in the slope at 10 18 35 eV (Unger et al. 
2007). Due to the uncertainties in the hadronic interac- 
tions at the highest energies, the interpretation of these 
elongation rate depends on the hadronic physics models 
used in the analysis, and is therefore rather ambiguous at 
present. However, regardless of which interaction models 
are used, the data appears to require the presence of a 
substantial fraction of heavy or intermediate-mass nuclei 
in the range of GZK cutoff energies. A possible conflict 
in the chemical composition between the above Auger 
elongation results and recent results by the Pierre Auger 
collaboration (2007) on large scale spatial correlations is 
an issue which remains to be resolved. 

The proposed astrophysical source models of UHECRs 
include AGN jets (e.g. Biermann 1987; Berezinsky et al. 
2006), intergalactic accretion shocks (e.g. Inoue, Aha- 
ronian & Sugiyama, 2005), and gamma-ray bursts (e.g. 
Waxman 1995, 2004a; Vietri 1995; Wick et al. 2004; Der- 
mer & Atoyan 2006; Murase et al. 2006). Recently, we 
have proposed that extra-galactic hypernovae associated 
with sub-energetic GRBs are also a probable source for 
UHECRs, whose energetics are sufficient large to account 
for UHECRs from the second knee and above (Wang et 
al. 2007). These UHECR accelerators also accelerate 
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electrons, which produce optical, x-ray or gamma-ray 
photons through synchrotron and inverse Compton emis- 
sion. Since a sufficiently high density of these low-energy 
photons can potentially lead to photo-disintegration of 
the cosmic ray heavy nuclei, a natural question is un- 
der what conditions do these sources allow the survival 
of UHE heavy nuclei from the acceleration site. In this 
paper we explore this problem for GRBs and hypernovae 
sources 5 , since they have the most compact source sizes 
and might therefore pose the greatest hurdle for UHE 
heavy nuclei survival. 

Long-duration GRBs are generally believed to result 
from the core collapse of massive stars, direct evidence 
for which comes from the spectroscopic identification of 
bright supernovae in association with these GRBs (see 
Woosley & Bloom 2006 for a review). The collapse leads 
to a black hole (or magnetar) plus accretion disk sys- 
tem. The annihilation of neutrinos or Poynting flows 
arising from the inner hot accretion disk produce a high 
entropy fireball outflow, which expands and converts its 
internal energy into the kinetic energy of a small amount 
of baryon material, while propagating through the star 
along the rotation axis of the collapsing core. The out- 
flow gets collimated and finally two highly relativistic jets 
break out of the stellar envelope. In the GRB scenario 
for UHECRs, both internal shocks (Waxman 1995) and 
external shocks (Vietri 1995; Wick et al. 2004; Dermer 
& Atoyan 2006) in the jets are suggested to be able to 
accelerate baryons to ultra-high energies. Internal shocks 
develop due to internal collisions between the shells in the 
fireball jet ejecta, while the external shocks occur when 
the fireball ejecta interacts with the surrounding inter- 
stellar medium (see Meszaros 2006 and Zhang 2007 for 
recent reviews). In the internal shock scenario, viewed 
as a source of the cosmic ray nuclei, the jets themselves 
must contain heavy nuclei, whereas in external shocks 
the heavy nuclei could come from the swept-up circum- 
burst material that the jet is running into, presumably 
the stellar wind or interstellar medium. Since the two 
shocks have quite different dissipation radii and emission 
properties, the disintegration problem of the UHE comic 
ray nuclei are different, as discussed in §£j4] and [5] respec- 
tively. 

Hypernovae are a peculiar type of supernovae with 
higher ejection velocities of the remnant shell and 
generally larger explosion energies than typical super- 
novae (Paczyhski 1998). The prototype of hypernovae 
is SN1998bw, a type Ic supernova associated with 
an under-energetic GRB, GRB980425 (Galama et 
al. 1998). Mildly relativistic ejecta components are 
inferred to be present in all three of the well-identified 
hypernova/under-energetic GRB systems so far, namely 
SN1998bw/GRB980425, SN20031w/GRB031203 and 
SN2006aj/GRB060218. Based on this mildly relativistic 
ejecta component and the event rates of these objects, 
we have shown that hypernovae can accelerate particles 
to ~ Z x 10 19 eV, and the energetics and occurrence rate 
is sufficient to account for the flux of UHECRs above 
the second knee (~6x 10 17 ) eV, where Z is the nuclear 
charge of the accelerated particles (Wang et al. 2007). 
The particles are accelerated in the hypernova blast 

5 The photo-disintegration of UHE nuclei in AGN scenario has 
been discussed in Anchordoqui et al. (2007b) and Dermer (2007a). 



wave formed by the interaction between the hypernova 
ejecta with the surrounding stellar wind medium. In §6] 
we study whether UHE heavy nuclei accelerated in such 
hypernovae can survive their photon environment. 

The goal of this paper is therefore to answer the ques- 
tion, if UHE heavy nuclei are accelerated in GRBs or 
hypernovae, what are the possible origins of these nu- 
clei, and under what conditions would they survive their 
source environment; or conversely, under what conditions 
would one expect these sources to accelerate mainly pro- 
tons to UHE energies. 

2. ORIGIN OF HEAVY OR INTERMEDIATE MASS NUCLEI 
IN THE SOURCES 

It is at present unknown whether heavy or intermediate 
mass nuclei are present in GRB jets. At the base of the 
outflow the jets start out as a hot fireball, within a region 
of size 10 6 — 10 7 cm and temperature kT = 1 — lOMeV, 
where any heavy nuclei will be photo-disintegrated due 
to the abundance of photons with energies comparable 
to the nuclear binding energy, ~ lOMeV. The fireball is 
thus initially made up of free nucleons, e ± pairs, trapped 
blackbody radiation and magnetic fields. 

As the fireball expands and cools, the free nucleons 
in the jet may recombine into a-particles (Beloborodov 
2003), but they will not form heavy nuclei. However, 
as the jet burrows through the stellar core (from in- 
side outward, Fe, O, C cores in sequence), heavy nuclei 
from the stellar surroundings could be entrained into the 
jet. According to numerical simulations of jet propaga- 
tion (Zhang et al. 2003), Kelvin-Hclmholtz instabilities 
and/or oblique shocks that develop lead to the mixing 
of surrounding material into the jet, while the jet is ad- 
vancing with a sub-relativistic velocity. Since the tem- 
perature of the thermal photons trapped in the jet de- 
creases significantly as jet expands, with T(r) ~ r , 
these thermal photons are no longer able to disintegrate 
the entrained nuclei when the jet has reached the radius 
of the stellar Fe core at ~ 10 9 cm. To avoid spallation 
of the nuclei due to nucleon-nucleon collisions, the rel- 
ative velocity between the jet and the surrounding core 
material should be below a critical value (3 sp ~ 0.14 (in 
units of the speed of light), at which the relative kinetic 
energy equals the nuclear binding energy ~ lOMeV. The 
calculations show that the jet head moves with a velocity 
about 10 9 cms _1 inside the Fe core at ~ 10 9 cm and then 
the velocity increases as ~ r 1 / 2 (Meszaros & Rees 2001). 
Thus, the Fe and O nuclei entrained from the surround- 
ings can survive both photo-disintegration and nuclear 
spallation. 

After the jet head breaks out of the star, in its subse- 
quent stages the jet continues to escape through the evac- 
uated funnel cavity through the star, so it accelerates to 
a relativistic velocity at radii still inside the star. There- 
after, the Kelvin-Helmholtz instability which causes the 
mixing of nuclei will be suppressed due to the relativistic 
relative motion (Zhang et al. 2003). Moreover, since the 
relative velocity between the jet and the surrounding star 
exceeds the critical velocity (3 sp , spallation effects will 
need to be taken into account. One finds that the time 
for spallation t sp — 1/ ! (a ap nNc) is much shorter than the 
dynamic time r/(Tc) when the relativistic jet is inside 
the progenitor star, where = L / (4irr 2 r 2 mNC 3 ) is the 
nucleon density of the jet and T is the Lorentz factor of 
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jet. This means that even if nuclei get entrained, they 
will be spallcd into lighter nuclei. Since the break-out 
time for initial jet after its inception is about % = 10 s, 
long GRBs with rest-frame durations t 1 ~ 2 — 10 s will 
have spent most of their life inside the star, while the jet 
head was moving sub-relativistically, so most of the jet 
can have become Fe-enriched before the internal shocks 
occur. For longer duration jets, with t 7 > 10s, only the 
first ~ 10 s portion of the burst may be Fe-rich, while the 
rest of the outburst may consist of pure protons. 

The above discussion is relevant to the models that in- 
voke the composition of the GRB jet as UHECR nuclei, 
which are accelerated in internal shocks that are widely 
assumed to be responsible for the prompt 7-ray emis- 
sion. For GRB external shocks and hypernova remnant 
blast wave shocks, the nuclei could be due to the ma- 
terial swept-up by the shock front. The progenitors of 
long-duration GRBs and hypernovae are thought to be 
Wolf-Rayet stars, as the spectral type of the discovered 
supernovae in these events is typically Ic. These stars 
are stripped of their hydrogen envelope and sometimes 
even the helium envelope. The heavy or intermediate 
mass UHE nuclei may originate from the stellar wind of 
the Wolf-Ryet star. In WC type Wolf-Rayet stars, the 
C abundance is Xq = 20% — 55% (by mass) and the O 
abundance is Xq = 5 — 10% (Crowther et al. 2007). In 
WO type Wolf-Rayet stars, O and C abundances are even 
higher, with X = 15- 25% and X c = 40% - 55% (e.g. 
Kingsburgh et al. 1995). The abundance of heavy ele- 
ments in stellar wind of these Wolf-Rayet stars is clearly 
much higher than the solar value. 

3. DISINTEGRATION OF UHE NUCLEI 

The most relevant processes that may prohibit acceler- 
ation of heavy nuclei to UHE at the astrophysical sources 
are the photodisintegration, photopion production and 
nuclear spallation. The former two processes are due to 
interactions of UHE nuclei with surrounding photons and 
the latter is due to interactions with other nuclei or nu- 
cleon. Both the photodisintegration and photopion pro- 
cesses are dominated by resonance production of cither 
an excited state of the nuclei or a Delta and subsequent 
de-excitation or decay. In case of photodisintegration 
this dominant channel is called a giant dipole resonance 
(GDR) in the energy range ~ 10 — 30 McV with a thresh- 
old energy of ~ 10 MeV per nucleon. In case of photopion 
the threshold energy for pion production is ~ 150 MeV. 

For an UHE nucleus with Lorentz factor ja propagat- 
ing through an isotropic photon background with differ- 
ential number density n{e) at energy s, the photodisin- 
tegration or photopion rate is given by (Stcckcr 1968) 
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where e' is the photon energy in the rest frame of the 
nucleus. The threshold energy eth = 10 MeV for pho- 
todisintegration and m^c 2 + m 2 c 2 /2m p ~ 145MeV for 
photopion. The cross-sections for photodisintegration in 
the energy range e t h < e' < 30 MeV with loss of one 
nucleon can be approximately described by a Lorentzian 
form (Puget, Stecker & Bredekamp 1976; Anchordoqui 
et al. 2007a) as 
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where Ag^jj and <jo,a are the width and maximum 
value of the cross section, so is the energy at which 
the cross section peaks. Fitted numerical values are 
(T , a = 1.45A x 10~ 27 cm 2 , A GDR = 8 MeV, and 
£ = 42.65.4~ - 21 MeV for A > 4 (Karakula & Tkaczyk 
1993). Above e' <~ 30 MeV and below the pion pro- 
duction threshold energy, photodisintegration may re- 
sult in multiple nucleon emission, although with much 
lower cross-section. Puget, Stecker & Bredekamp (1976) 
suggested a parametrization such that the cross-section 
(TA(e') integrated in the range Eth < e' ;$ 30 MeV and in 
the range 30 MeV < e' < 150 MeV are equal. We assume 
a flat crosss-section in the range 30 McV < e' < 150 MeV 
of 6.6 x 10~ 27 cm 2 and 1.7 x 10~ 27 cm 2 respectively for 
iron and oxygen nuclei, satisfying above condition. This 
more acccurate cross-section affects the photodisintegra- 
tion rate in case the photon spectrum is very hard as 
we will see later. For soft photon spectrum equation 
(2) is adequate and often a delta function approxima- 
tion <ja(£') ~ 0o, A^GDRS(e' — so) may provide order of 
magnitude estimate. 

The photopion cross-section formula for Delta reso- 
nance production is well known and is parametrized by 
Miickc et al. (2000) as 
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where s = m 2 c 4 + 2e'm p c 2 is the center-of-mass energy, 
00, 7T = 3.11 x 10~ 29 cm 2 and the peak cross-section is 
- 4.12 x 10~ 28 cm 2 at s' w 0.3 GcV. The width of the 
resonance is A T = 0.11 GeV. The photonuclear cross 
section is a v ~ A 2 ^a/^. Comparing this cross section 
to that of photo-disintegration for an iron nucleus with 
A = 56 or Oxygen nucleus with A = 16, we can see that 
the photopion cross section is smaller than the photo- 
disintegration cross section. As a result, the photodis- 
integration rate is typically higher than the photopion 
rate. 

The rate for spallation process due to collision with 
other nucleons in the sources may be estimated as 
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where cr sp = 5 x 10- 26 A 2 ^ cm 2 is the cross section for 
spallation of a nucleus with atomic number A, n?q is the 
number density of nucleus in the sources and (3 is the 
velocity of the UHE nucleons. As we show below, even 
for internal shocks which have the highest nucleon den- 
sity, this spallation effect is less important than photo- 
disintegration. 

Below we will discuss the photodisintegration process 
in all three possible UHE nuclei sources, i.e. GRB inter- 
nal shocks, GRB external shocks and hypernova remnant 
blast waves, respectively. 

4. GRB INTERNAL SHOCKS 

Due to the variable injection at the base, internal col- 
lisions occurr within the unsteady plasma jet which de- 
velop into internal shocks. The prompt gamma-ray burst 
emission is supposed to result from the non-thermal emis- 
sion of electrons accelerated in these shocks. The two 
main processes through which the heavy nuclei could be 
disintegrated in internal shocks are photo-disintegration 
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or photopion interaction by the X-ray photons in prompt 
emission, and spallation due to collisions with other nu- 
clcons in the rclativistic jet. 

4.1. Nucleus photo- disintegration due to prompt X-rays 

The GRB prompt photon spectrum is well fitted in 
the BATSE range (10 KeV - 3 MeV) by a combina- 
tion of two power-laws, n(e) oc e _/3 , with different val- 
ues of (3 at low and high energy (Band et al. 1993). 
The break energy £b,obs in the observer frame is typically 
£b,obs = Teb ~ lMcV, with (3 ~ 1 at energies below the 
break and (3 ~ 2 above the break, where T is the bulk 
Lorentz factor of the relativistic flow that develops inter- 
nal shocks and e& is the break photon energy in the flow 
rest-frame. At low energy, in the framework of the inter- 
nal shock model, there could be another spectral break, 
i.e., the synchrotron self- absorption break e ssa . 

Thus the photon spectrum, in the comoving frame, is 
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where rib = n(eb) is the photon number density at the 
break energy in the comoving frame of the wind. 

Using this photon spectrum, the inner integral in 
Eq.(l) gives j~ 2i n{e)e-*de = ^(gi-J-d+fl. The 
energy density of photons in the wind rest-frame in the 
energy range of BATSE is U-y = meKl /2 + ln(e m / £b)) ~ 
2nbS 2 , where sm — 3MeV is the upper energy range 
of BATSE. This energy density is related to the ob- 
served photon luminosity of GRBs, whose typical value 
is L 7 ~ 10 51 ergs _1 , by L 7 = AiTR 2 n T 2 cU 1 , where R in is 
the radius of the emitting region of these photons. 

The first integral of Eq.(l) can be done numerically for 
the exact cross sections in Eq.(2) nd (3). For the sake of 
an analytic treatment, we approximate the cross section 
as being mainly contributed by the resonance peak and 
find that 
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where k = 2je SS a/so- The observer- frame energy of the 
nucleus that interacts preferentially with photons with 
energy Sb and e ssa are, respectively 
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where T ~ 300r 2 .5 is a typical GRB bulk Lorentz factor. 

Now we can obtain the time tdis of an nucleus of energy 
E in the observer frame: for Eb < E < E ssa 
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The dynamic time of internal shock is 

^dyn = Rin/cT — lR in ^T 2 \s. 



(8) 
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By comparing Eq.(8) or Eq.(9) with Eq.(10), we can see 
that the heavy nuclei can survive when either 
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in the case of idis S tdyn- 

For comparison, the acceleration time of UHE nuclei 
with nuclear charge Z in the internal shock is 
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where a ~ few, describing the ratio between the ac- 
celeration time and Larmor time, B = ( j^sfj^- ) 1 / 2 = 
8 x lQ 4 e^ 2 _ l L 1 J^ 2 R^Y 2 \G is the comoving frame mag- 
netic field in internal shocks. The synchrotron loss time 
for UHE nuclei is 

67rm^c 3 r A , 

Wn= a T m 2 EB 2i Z 

(14) 

which is a factor of (A/Z)* ~ 16 longer than that of 
protons. 

By equating i acc with idyn or i syn , we get the maxi- 
mum energies of accelerated nuclei for these two cases 
respectively, 
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We compare these timescales in Fig.l for two differ- 
ent sets of parameters, i.e. R = 10 13 cm, T — 10 2 ' 5 (top 
panel) and R = 10 14 cm, T = 10 3 (bottom panel). We can 
see that: 1) the iron nucleus can be accelerated to ener- 
gies > 10 20 eV, since in both cases, i acc < £dyn < isyn for 
E = 10 20 eV; 2)when the internal shock radius and/or 
the bulk Lorentz factor are relatively large, the photo- 
disintegration time of iron nuclei is longer than the dy- 
namic time, as seen in the bottom panel of Fig.l. Note 
that the photo-disintegration time and the photopion in- 
teraction times in these figures are calculated numerically 
according to Eq.(l) with exact cross sections. 

For smaller internal shock radii (Fig.l, top panel), the 
photo-disintegration effect becomes important if the pho- 
ton spectrum n{e) ~ e -1 extends to < IKeV. This is 
consistent with the result in Anchordoqui et al. (2007b). 
However, if there is a self-absorption break at several 
KeV in the photon spectrum that leads to a drop in 
the number density of the photons with which the UHE 
nucleus mainly interacts, the photo-disintegration effect 
will be suppressed, as can be seen from the blue dashed 
lines, where a synchrotron self-absorption break at 5 KeV 
in the photon spectrum has been assumed. This is pos- 
sible for a small internal shock radius, since a smaller in- 
ternal shock radius might lead to higher synchrotron self- 
absorption break e ssa ■ The value of e ssa also depends on 
the dissipation model as well as the radiation mechanism 
for prompt emission (e.g. Panaitescu & Meszaros 2000; 
Rees & Meszaros 2005; Pe'er et al. 2005; Pe'er & Zhang 
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2006), which is, however, largely uncertain at present. 
For the simple optically thin internal shock model, e SS a 
is estimated to be 



-ssa,obs 



(16) 



For UHE intermediate-mass nuclei, such as O nuclei, 
the constraints for survival, Eqs.(ll) and (12), are much 
looser. Thus, for small internal dissipation or shock radii, 
intermediate-mass nuclei are much more likely to survive 
photo-disintegration. This can be also seen in Fig. 2, 
where all the timescales are calculated for O nuclei with 
Z = 8 and A = 16. 

4.2. Nuclear spallation due to collisions 

An UHE nucleus can also spallate due to collision with 
other nucleons in internal shocks. The nucleon den- 
sity of a relativistic flow with kinetic energy luminosity 
Lk = 10 52 ergs _1 at internal shock radius Ri n is njy = 
L k /(^R 2 n T 2 m N c 3 ) = 1.8 x 10 12 L fc , 52 i^ 13 r^ciii- 3 . 
So the spallation time is 

1 



t 



sp 



nNO- S pC 



r 2 (^_)-2/3 
fe,52-"'m,13 i 2.5lgg7 b ' 



25L7l R 2 



(17) 



Comparing this with the dynamic time, we find that the 
heavy nucleus can survive if 
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This condition is much easier to satisfy, compared to 
the more stringent photo-disintegration constraint of 
Eq.(ll). 

5. GRB EXTERNAL SHOCKS 

As the relativistic GRB jets expand and sweep up the 
surrounding material, external shocks develop. The ac- 
celerated electrons in the external shock produce the ob- 
served afterglow emission. The external shock has also 
been suggested to be able to account for UHECRs. The 
early deceleration phase, when the jet converts about 
half of its energy to the swept-up material, is the most 
effective phase for the external shock to accelerate par- 
ticles, since the shocked material has the largest Lorentz 
factor at that time. The early external shock produces 
the early x-ray afterglow, whose average luminosity ob- 
served by Swift is about L x ~ 10 48 ergs _1 during the 
first ~ 10 - 100 s (Nousek et al. 2006). In this external 
shock scenario, the accelerated UHE nuclei with energy 
E > 10 19 eV preferentially interact with these x-ray pho- 
tons, if the Lorentz factor is typically 10 2 < T < 10 3 . 

Taking the early X-ray afterglow spectrum in the fast- 
cooling regime (e.g. Meszaros 2006) with F v ~ 
corresponding to photon spectral index (3 — 2, we can 
obtain the photo-disintegration rate of a nucleus moving 
with Lorentz factor 7^4, i.e. 
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where U x is the comoving-frame energy density of X-ray 
afterglow photons and k = ln(ex,M/ex,m) — 3, ex,M and 
cx,m being the upper and lower end of Swift XRT energy 
threshold. The energy density is related to the luminosity 
by L x = 4:TtRI x T 2 cU x , where R ex is the radius of the 
external shock at the end of the free expansion phase of 



the ejecta. Now we can obtain the photo-disintegration 
energy loss time of a cosmic ray nucleus of energy E 
propagating in the early afterglow photons 

t dis = 3 x io 6 i4 il7 ri. 5 L- 1 48 ^ 1 (^)-°- 21 s. (20) 

Comparing this with the dynamic time in the comov- 
ing frame of external shocks td yn = R e x/(2cT) = 
10 4 i? ea; ,i7r^5S, we can see that the heavy nucleus can 
survive when 

Re X ,nTi 5 > 3 x lO-^^o^) - 21 . (21) 

For a constant density medium environment , this gives 
an upper limit on the number density of the medium 

n < 2.5 x 10 6 S fc ,53r^ 5 L- 3 48 S 20 3 (A/56)- - 63 cm- 3 , (22) 

while for a stellar wind environment with a density pro- 
file p = (4^-) r — 2 = 5 x 10 1:L A*r _2 gcm _1 , this gives 
constrains on the mass-loss rate parameter 

A, < l.ZE^r^L-^E^iA/W)- - 21 , (23) 

where A+ = 1 corresponds to M = 10" 5 M Q yr _1 and 
v w = lOOOkms^ 1 . 

After the initial free expansion phase, the external 
shock starts to decelerate as more and more material 
is swept-up. The shock radius evolves with time as 

Rex{t) = 4 x 10 17 E^^Uq 1 ' a t 1 ^ 4 cm for constant density 

medium and R ex (t) = 2 x 10 18 £^ ^A'^t^^cm for 
wind medium, where td is time in units of days. Assum- 
ing that the magnetic field energy density acquires a frac- 
tion 6b — 0.1 of the internal energy, the magnetic field 
evolves as B 



-, 1/2 ,-,1/8 3/8,-3/8^, r 

l^B -i-kfc 53 n o t d G for a constant 



density medium and B = 0.07e^' 2 _ 1 £' A . ^ 3 4 A^^ 1 t d 3 ^ 4 G 
for wind medium. Since the synchrotron loss is typically 
unimportant for UHE nuclei in external shocks, the max- 
imum energy of particles accelerated in external shocks 
is obtained by equating i acc with tdyn, which gives 



£max — ZeBR ex (t) 

-121/2.^1/2 p 3/8 l/ u 
J \26) t B.-l rj k.53 n L d 
1Q 21 ( Z )e l/2 E l/4 A l/4 ,-1/4 



3 x 10 2 



* 26 ^ c B,-l^fc,53 



it. 



eV (ISM) (24) 
eV (wind) 



This shows that the maximum energy of accelerated 
particles by external shock decreases rather slowly with 
time, so that even weeks to months after the burst, heavy 
nuclei can be still accelerated to UHE energies by after- 
glow shocks, even if the flux may decrease with time. 

The optical depth r for photo-disintegration of these 
nuclei evolves with time as 



tdyn 
^dis 



L X T 4 R ex \ t 



t 1 / 4 (ISM) 
i-V 2 (wind) 



(25) 



where an X-ray afterglow luminosity L x ~ t^ 1 has been 
assumed. The slowly increasing optical depth or decreas- 
ing optical depth indicates that the accelerated heavy nu- 
clei by GRB external shocks can survive for a relatively 
long time. We plot in Fig. 3 the evolution of photo- 
disintegration optical depth and the maximum energy 
of accelerated particles for iron nucleus (top panel) and 
Oxygen nucleus (bottom panel) respectively. 
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6. HYPERNOVA REMNANT BLAST WAVE 

The observations of the radio afterglow of the hyper- 
nova SN1998bw showed that about 10 50 erg of kinetic 
energy were released in the form of a mildly relativis- 
tic ejecta (Kulkarni, et al. 1998; Chevalier & Li 1999). 
The interpretation of the X-ray afterglow also favors a 
mildly- relativistic ejecta component (Waxman 2004b) . A 
recently detected strong thermal X-ray emission compo- 
nent in another sub-energetic burst (GRB060218), as- 
sociated with SN2006aj, may also be associated with a 
mildly- relativistic supernova shock breakout, in which 
the mildly relativistic supernova ejecta has an energy 
> 10 49 erg (Campana et al. 2006). Due to the large su- 
pernova explosion energy and the much lower than typ- 
ical GRB energy, attempts have been made to ascribe 
the prompt gamma-ray emission to the shock from the 
mildly relativistic ejecta as it breaks out through the hy- 
pernova progenitor's outer envelope (Woosley et al. 1999; 
Tan et al. 2001), although a generally accepted conclu- 
sion has not yet been reached. We have used the term 
semi-relativistic hypernovae to denote such supernovae 
exhibiting a mildly relativistic ejecta component, seen 
in association with GRBs. Such high-velocity ejecta are 
a key ingredient of the hypernova model for UHECRs, 
since the low-velocity bulk ejecta is not able to acceler- 
ate particles to such high energies. Based on the shock 
breakout scenario, we have suggested that there might 
be a continuous distribution of ejecta energy in velocity, 
i.e. E k (T(3) ~ 3 x 10 52 (r/?/0.1)- Q with a ~ 2, where (3 
is the ejecta velocity in units of speed of light (Wang et 
al 2007). 

As the hypernova ejecta expand, they transfer their 
energy to the swept-up stellar wind, and external shocks 
develop, which accelerate wind particles to ultra-high en- 
ergies. Since the ejecta has a velocity distribution pro- 
file, the leading edge higher-velocity ejecta decelerate the 
earliest and then the lower-velocity ejecta decelerate pro- 
gressively. The maximum energy of accelerated particles 
is related to the ejecta velocity T(3 by 



~ ZeBRf3 



= 1.3 x lO^f)^ ($) Al^eV. 



(26) 



Although higher velocity ejecta can accelerate particles 
to higher energies, the kinetic energy in ejecta of Y[3 > 
0.5 are too low to account for the UHECR flux (Wang et 
al. 2007), due to the steep distribution of energy E k oc 
(r/3)~ 2 . For ejecta with a velocity F/3, the free expansion 
phase before deceleration sets in lasts for a time 

t f = 1300(r/3/0.5)" 5 A- 1 days (27) 

and the radius of the ejecta at this time is 

r f = 1.7 x 10 18 (r/3/0.5)" 4 A- 1 cm. (28) 

There are two photon sources which could cause photo- 
disintegration of heavy nuclei: one is provided by hyper- 
nova thermal photons from radioactive elements of the 
hypernova ejecta, and another is the synchrotron pho- 
tons from the hypernova remnant blast wave. For the 
first mechanism, we use the luminosity of SN1998bw as 
a representative. At time t ~ 1000 days after the burst, 
the optical luminosity of SN1998bw drops to the level of 
about L H e ~ lO^crgs- 1 (Sollcrman et al. 2002). A nu- 
cleus of energy E — 10 20 eV interacts with target photons 



with energy e t > 0.01(A/56)£; 2 - 1 eV. A rough estimate of 
the optical depth of photo-disintegration of heavy nuclei 
due to hypernova thermal photons is 



< 



„ ( ^HE 

u V 4?rr^c£HE 



3x lO-^HE^/jsdgf)- 1 , 



(29) 



where r\ ~ 4 is the compression ratio of the hypernova ex- 
ternal shock and £he — leV is the characteristic energy 
of hypernova thermal photons. We can see that thermal 
photons are so sparse that they have a negligible effect 
on the photo-disintegration of UHE nuclei. 

The synchrotron emission from the stellar wind 
shocked by the hypernova remnant could be brighter 
than the hypernova emission itself at the late stages of 
the hypernova, e.g. tf ~ 1300days after the burst in 
our case, since the hypernova luminosity drops rather 
quickly after the peak. We estimate the luminosity 
from the shocked wind at the time when the T[3 = 0.5 
ejecta begins to decelerate. The total number of the 

swept-up wind particles is N e — (^r") r / = 6.6 x 

lO^^/lSOOd)" 4 / 5 ^ 475 , where d denotes time in days. 
Assuming that the magnetic field energy density is am- 
plified to a fraction es of the shock internal energy, 
the magnetic field is B = 0.015(r/3/0.5) 5 A^ /2 G = 
O-O^iz/lSOOd)- 1 ^ 7 ^. Assuming a power-law en- 
ergy distribution for accelerated electrons dn e /dj e cx 7~ p 
with p ~ 2, one obtains the characteristic frequency 
v m = 7 x 10 5 e 2 _ 1 (i / /1300d)- 9 / 5 A* 3/10 Hz, where e e is 
the fraction of the shock energy that goes into electrons. 
High energy electrons will cool in the magnetic field and 
cause a break in the electron energy distribution. The 
characteristic synchrotron frequency corresponding to 

this break is v c = 4 x 10 13 (t/ /1300d)yL7 3/2 Hz. The lumi- 
nosity at the frequency v m is L Vm — N e ^ m P(^ e ^ m ) ~4x 
10 36 e 2 _ 1 e B ,_i(t / /1300d)- 2 A*crgs- 1 . Forp~2, the lu- 
minosity at frequency v is L y — L v 

{v/v m ) 1/2 for v m < 
v < u c , while for v > v c , L v ~ L Vc — L Vm (y c jv m ) l l' 1 . 
The luminosity at energy e t — 0.01(A/56)E^ ) 1 eV, cor- 
responding to the energy of the photons with which the 
UHE nuclei of energy E interact at the resonance peak, is 
L syn = 10 40 e e ,_ie 3 3 /4 _ 1 (i / /1300d)- 11 Ai 15 ergs- 1 . So we 
see that for wind parameters A+ = 1, the synchrotron 
luminosity does exceed the hypernova luminosity at the 
time ~ 10 3 days 6 . Once we know the synchrotron lumi- 
nosity of a hypernova remnant and the photon spectral 
index {[3 — 3/2), we can, using Eq.(l), get the optical 
depth of photo-disintegration of heavy nuclei due to such 
synchrotron photons, i.e. 

T ^o vlft -3, 3/4 ,1.35 W l/2,A ^ n tf 19 



Using t < 1, we find that, after a time 



(30) 



A 



tf > l00e™\e%U^E% 26 (g)°- 37 d, (31) 

6 Unlike in SN2006aj, the inferred stellar wind for SN1998bw is 
much weaker, with A* = 0.04 — 0.1, which may explain why the 
optical emission of SN1998bw, detected at ~ lOOOdays after the 
burst, is still contributed by the hypernova emission 
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the nucleus can survive in the hypernova synchrotron 
photon environment. Note that the free expansion time 
for ejecta with T/3 < 0.5 (Eq.27) is longer than this time. 
The maximum energy of accelerated particles, Eq.(26), 
depends on the time tf as 

e max ~ 1.3 x 10 20 ( -) e 1 ' 2 . ( -^-] A^eV. 
\26J B ~ 1 V1300dy 

(32) 

The evolution of the optical depth for photo- 
disintegration and the maximum energy of accelerated 
nuclei are shown in Fig. 4. 

7. DISCUSSION AND CONCLUSIONS 

The possible presence of heavy nuclei in UHECRs 
raises interesting questions: what is the origin of these 
nuclei? can these nuclei survive in the sources where 
they get accelerated? In this paper, we have endeav- 
ored to address these questions for two proposed UHECR 
sources discussed in the literature, namely GRBs and hy- 
pernovae. For GRBs, both internal shocks and external 
shocks have been suggested to be able to accelerate par- 
ticles to UHE energies, and in this paper we have consid- 
ered the role of both of these in the context of UHECR 
heavy nuclei. We have sketched out some possible mech- 
anisms for the presence of heavy nuclei in these scenarios. 
For GRB internal shocks, we suggest that the nuclei are 
entrained from the progenitor stellar core, e.g. the Fe 
core, the O core, etc., during the stage when the jets 
arising from accretion in the innermost collapsing core, 
are making their way out through the star. Since in- 
stabilities favoring cntrainment are predominant mainly 
during the initial stellar crossing phase of the jet, a jet 
composition with a substantial fraction of heavy nuclei 
may only be expected in long bursts with observed 7-ray 
durations t 7 < 10 — 15(1 + z) s, while longer burst jets 
would be expected to consist mainly of protons. For the 
GRB external shock and hypernova source scenarios, the 
nuclei may be the heavy elements present in the stellar 
wind of Wolf-Rayet stars, which arc thought to be the 
progenitors of GRBs and hypernovae. The stellar wind 
of Wolf-Rayet stars, especially the WO, WC sub- types 
stars, is heavily enriched with intermediate mass nuclei, 
such as O, C, etc. 

After they escape from their sources (see e. g. Der- 
mcr 2007b), UHECR nuclei can be subject to photo- 
disintegration also in intergalactic space, before arriv- 
ing at Earth. An UHE heavy nucleus with an energy 
> 10 19 eV will mainly collide with cosmic infrared back- 
ground (CIB) photons. Using new constraints on CIB 
data, it has been found that UHE iron nuclei with en- 
ergy < 10 20 eV have a mean free path of > 300Mpc and 
that this mean free path increases as the energy of the 
nucleus decreases (Hooper et al. 2007; Stecker & Sala- 
mon 1999). Lighter nuclei have relatively shorter mean 
free path (Hooper et al. 2007). Thus, UHE nuclei can in 
principle originate from sources at cosmological distance, 
such as GRBs and hypernovae. The main question ap- 
pears to be whether they survive the environment of their 
original sources. For this reason, most of the present pa- 
per has been devoted to a quantitative study of the sur- 
vival of UHE nuclei in the above mentioned sources. We 
find that: 

i) In GRB internal shocks, heavy nuclei can survive 



in the sources if the internal shock radius and/or the 
Lorentz factor of the relativistic jets are relatively large, 
as given by Eq.(ll). Thus one might expect acceleration 
of Fe nuclei to be more favored in bursts with smoother, 
longer variability timescale light-curves. For a smaller in- 
ternal shock radius, the photodisintegration process due 
to prompt X-ray photons may become optically thick. 
However, if a synchrotron self-absorption break is present 
above several keV in the photon spectrum, the reduced 
number density of X-ray photons will lower the photodis- 
integration optical depth accordingly. In general, com- 
pared to heavy nuclei, the UHE intermediate-mass nuclei 
find it easier to survive photo-disintegration in internal 
(and external) shocks. 

ii) In GRB external shocks, due to the much larger 
dissipation radii compared to internal shocks, UHE nu- 
clei can easily survive in the sources. We have also cal- 
culated the evolution of the photo-disintegration optical 
depth during the afterglow phase for both an ISM ex- 
ternal medium and a stellar wind medium. The results 
show that UHE nuclei can survive in the afterglow shock 
for a relatively long time in both cases. 

iii) In the hypernova remnant acceleration scenario, 
UHE nuclei can survive in the sources, except in the early 
short period of time (< lOOdays for typical parameters, 
see Eq.(31)). In this early short period, however, only a 
very small amount of ejecta energy has been converted 
to UHECRs, so their contribution to the UHECR flux 
is negligible. Most of the UHECRs originate from the 
blast wave ejecta with T(3 < 0.5, which decelerate after 
~ lOOOdays, typically. UHE nuclei accelerated during 
this time are safe from photo-disintegration. 

Earlier calculations on the photo-disintegration of 
UHE nuclei in GRB internal shocks by Anchordoqui et 
al. (2007b) consider only one specific set of typical pa- 
rameters for the internal shock radius and the relativistic 
Lorentz factor. We have improved this by considering the 
whole parameter space and working out the constraints 
on the physical conditions under which UHE heavy nuclei 
can survive in internal shocks. We have also considered 
the effect of a self-absorption break in the prompt emis- 
sion on the photo-disintegration problem. Besides in- 
ternal shocks, we have explored the photo-disintegration 
problem of UHE nuclei in GRB external shocks and hy- 
pernova remnants. 

In summary, we have suggested possible scenarios 
for the injection of heavy nuclei into the acceleration 
zones of extragalactic sources such as GRB and hyper- 
novae. We found that, if heavy nuclei are accelerated 
in these sources, they will survive the threat of photo- 
disintegration under fairly general conditions for the case 
of GRB external shocks, and for hypernovae. They could 
survive also in GRB internal shocks, if the latter occur at 
relatively large radii and/or the bulk Lorentz factors are 
large. On the other hand, for small shock radii and/or 
smaller bulk Lorentz factors, a pure proton UHECR com- 
position would be favored if the self-absorption break in 
the photon spectrum is not high. Since the instability- 
induced entrainment process of heavy nuclei in internal 
shocks is currently not well-known, we stress that the 
expected fraction of heavy nuclei injected into the inter- 
nal shock acceleration process is uncertain. A significant 
fraction is plausible, but if the entrainment is inefficient, 
we would in any case expect a proton-dominated compo- 
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sition in GRB internal shock models. In GRB external 
shocks and hypernova shocks, on the other hand, the 
abundance of heavy nuclei is dependent on the external 
medium or the stellar outer envelope and wind enrich- 
ment fraction. 

As we were completing our work, the Pierre Auger 
Collaboration (2007) reported a plausible correlation be- 
tween UHECR at energies above 6 x 10 19 eV, assum- 
ing that they are protons, and AGNs within 75 Mpc se- 
lected from a particular catalog. Although this result 
is statistically significant at the 2.8a level, as the au- 
thors themselves stress, it does not rule out a possible 
origin in sources, e.g. galaxies, which are distributed 
in a similar manner as these AGNs. There is an unre- 
solved tension between the above Auger spatial corre- 
lation analysis suggesting protons, and previous Auger 
results on maximum shower elongations X max suggest- 
ing a significant heavy element component for UHECR 
in the same energy range. The AGN hypothesis for 
the UHECRs has been questioned in the analysis of Gor- 
bunov et al. (2007). In our present paper, we have inves- 
tigated the conditions under which UHECR originating 
from GRB or hypernovae could contain a significant frac- 
tion of heavy nuclei, as well as the conditions under which 
these UHECR would be expected to be mainly protons. 
If any of the observed UHECR were heavy nuclei, these 
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Fig. 1. — Top panel: Timescales in the internal shock (or in general, small dissipation radius) scenario, comparing the dis- 
integration time idis, the photopion interaction time t v , the dynamic time tdyn, the particle acceleration time t a cc and the 
synchrotron loss time t syn for iron nuclei at different energy E. The disintegration time (thick blue lines) and photopion interac- 
tion (thick red lines) times are calculated numerically by using exact cross-sections, i.e., equation (2) and multi-nucleon emission 
cross-section for photodisintegration and equation (3) for photopion. The thin blue lines corresponds to photodisintegration 
time using equation (2) only, i.e., ignoring the multi-nucleon emission cross-section. For the photodisintegration and photopion 
times, the solid line and dashed lines are for a synchrotron self-absorption break e 3sa ,obs at 1 KeV and 5 KeV respectively. The 
dissipation radius and jet Lorentz factor used in this plot are Ri n = 10 13 cm and F = 10 2 ' 5 . Note that the maximum energy of 
an accelerated iron nucleus is ~ 3 x 10 21 eV from t acc = t Byn at this radius. Bottom panel: The same as above but for a larger 
dissipation radius, R\ n = 10 14 cm and a larger jet Lorentz factor F — 10 3 . The maximum energy of an accelerated iron nucleus 
is ~ 2 x 10 21 eV from i acc = tdyn- The acceleration time has been roughly approximated as the Larmor time, i.e. a — 1 in both 
plots. 
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Fig. 3. — External shocks scenario. In the top panel, dashed lines show the evolution of the maximum energy of an accelerated iron 
nucleus for an external shock moving into a constant density medium (thick dashed line) and into a wind medium (thin dashed line) , while 
the solid lines show the evolution of the photo-disintegration optical depth of an iron nucleus of energy E = 10 20 eV in external shocks 
moving into the constant density medium (thick solid line) and the wind medium (thin solid line). The parameters used for the constant 



density medium case and the wind medium case are, respectively, {n = 1cm 



10^ erg, T 



10 x 



E = 10 53 crg, T = 10 2 - 5 , e B = 0.1 }. Bottom panel: the same as the top panel, but for UHE Oxygen nuclei. 



e B = 0.1 } and {A* = 0.1, 
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Fig. 4. — Hypernova acceleration scenario. Dashed lines: evolution of the maximum energy of an accelerated iron nucleus (thick dashed 
line) and an Oxygen nucleus (thin dashed line) in the hypernova remnant; Solid line: evolution of the photo-disintegration optical depth of 
an iron nucleus (thick solid line) and an Oxygen nucleus (thin solid line) of energy E = 10 20 eV. The parameters used are A* = 1, e e = 0.1 
and £b =0.1. 



